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ABSTRACT 

Context. Complex bipolar shapes can be generated either as a planetary nebula or a symbiotic system. The origin of the material 
ionised by the white dwarf is very different in these two scenarios, and it complicates the understanding of the morphologies of 
planetary nebulae. 

Aims. The physical properties, structure, and dynamics of the bipolar nebulae, M 2-9, Mz 3, and Hen 2-104, are investigated in 
detail with the aim of understanding their nature, shaping mechanisms, and evolutionary history. Both a morpho-kinematic study and 
a spectroscopic analysis, can be used to more accurately determine the kinematics and nature of each nebula. 

Methods. Long-slit optical echelle spectra are used to investigate the morpho-kinematics of M 2-9, Mz 3, and Hen 2-104. The 
morpho-kinematic modelling software SHAPE is used to constrain both the morphology and kinematics of each nebula by means of 
detailed 3-D models. Near-infrared (NIR) data, as well as optical, spectra are used to separate Galactic symbiotic-type nebulae from 
genuine planetary nebulae by means of a 2MASS J-H/H-Ks diagram and a AA'ih'ifti.y vs. TSOOT/HjS diagnostic diagram, respectively. 
Results. The best-fitted 3-D models for M 2-9, Mz 3, and Hen 2-104 provide invaluable kinematical information on the expansion 
velocity of its nebular components by means of synthetic spectra. The observed spectra match up very well with the synthetic spectra 
for each model, thus showing that each model is tightly constrained both morphologically and kinematically. Kinematical ages of 
the different structures of M 2-9 and Mz 3 have also been determined. Both diagnostic diagrams show M 2-9 and Hen 2-104 to fall 
well within the category of having a symbiotic source, whereas Mz 3 borders the region of symbiotic and young planetary nebulae in 
the optical diagram but is located firmly in the symbiotic region of the NIR colour-colour diagram. The optical diagnostic diagram 
is shown to successfully separate the two types of nebulae, however, the NIR colour-colour diagram is not as accurate in separating 
these objects. 

Conclusions. The morphology, kinematics, and evolutionary history of M 2-9, Mz 3, and Hen 2-104 are better understood using the 
interactive 3-D modelling tool shape. The expansion velocities of the components for each nebula are better constrained and fitted 
with a vector field to reveal their direction of motion. The optical and NIR diagnostic diagrams used are important techniques for 
separating Galactic symbiotic-type nebulae from genuine planetary nebulae. 

Key words. Planetary nebulae: individual: (M 2-9, Mz 3, Hen 2-104) - Stars: binaries: symbiotic - Stars: kinematics and dynamics 
- Stars: winds, outflows - Stars: jets - Infrared: stars 


1. Introduction 


The formation of complex bipolar optical nebulae, such as M 2- 
9, Mz 3, and Hen 2-104, continue to intrigue and perplex as¬ 
tronomers. There are competing theories to explain the bipolar 


morphologies observed: common envelope evolution (Nordhaus 
& Blackman||2006 1 or intrinsically bipolar ejections (Balick & 


Frank|2002|l. Other features, such as knots or clumps, are most 


often associated with interactions between slow-moving ejecta 
and a later fast wind (|Sadakane et al.|2010|l. 


A planetary nebula can form when a low mass star evolves 
through the AGB phase and subsequent thermal pulses, exposing 
a C-O white dwarf (WD) core with an effective temperature high 
enough that Lyman continuum photons ionise the ejected en¬ 
velope. Planetary nebulae (PNe) are observed to exhibit a wide 
variety of morphologies and in general a binary system is not 
a prerequiste to forming a PNe. However, for strongly bipolar 
planetary nebulae (bPNe), often with collimated outflows and 
point symmetric structures a binary central system is widely be¬ 


lieved to be essential for the extreme deviation from spherical 
symmetry seen. 

In a binary system of two low mass stars, when the orig¬ 
inal PN from the higher mass primary star has dispersed, the 
secondary star may then in turn evolve away from the main se¬ 
quence and form a second planetary nebulae. For widely spaced 
binary stars, the second PN may appear a little different to that 
from a single star. A symbiotic binary system (hereafter SySt), 
however, is a long-period interacting binary system composed 
generally of a hot white dwarf and an evolved mass-losing cold 
red-giant star, usually a Mira variable. In SySts, mass that is 
transferred and ejected from the secondary can be ionised by the 
old primary WD, forming another, different kind of PN termed 
a symbiotic nebula that shares many of the overall features of a 
classical PN. 

PNe and SySts are at distinctly different stages of stellar evo¬ 
lution. Since a SySt consists of a mass-losing Mira (or a red giant 
that has not yet reached the PN stage) and a white dwarf compo¬ 
nent that has evolved past the PN stage, they can be described as 
being simultaneously post-PN and pre-PN (Lopez et al.||2004}. 
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While SySts are binaries by definition, PNe may occur from ei¬ 
ther single- or multiple-star systems. It is not obvious whether 
a given bipolar nebula is a classical planetary nebula or a sym¬ 
biotic nebula. Symbiotic nebulae are sometimes misclassified as 
PNe, since both types of objects share the same characteristics 
(i.e., morphology and ionization structure). Additional evidence 
is required to identify a SySt. 

The optical spectra from symbiotic nebulae show both ab¬ 
sorption and emission lines from the surface of the cool red 
giant and the surrounding nebula, respectively. This requires 
the presence of a hot compact source of ionising radiation. 
Since both PNe and symbiotics differ vastly in their density, the 
change in intensities of their [O iii] emission lines at 4363 A and 
5007 A can provide a means of separating them. O vi Raman 
scattered lines at 6830 A, 7088 A seen only in symbiotic nebulae 
due to the extreme high densities in their cores ( Schmid|1989 1 . 
Broad-winged, single, or double-peak Ha line profiles, also indi¬ 
cate the presence of a symbiotic core ( Van Winckel et al. [1993) 1. 
It can also be possible to directly detect the contribution of a cool 
giant companion in the red part of the optical regime. The con- 
tinua are characterised by a dominant red continuum, typical of 
an F-M type giant for SySts, and a weak atomic recombination 
continuum in the optical for PNe. 


The NIR colour-colour diagrams (Whitelock & Munari 


1992[[Sdimeja & Kimeswenger|200H|Phillips|2007HSantander- 


et al.|2001t|Crocker et al.|200^|Corradi|2003t[Sokoroski|20031 


Corradi et al.|201 1[ [Werner et al j2014|l. 

To understand the origin of the range of morphologies of 
PNe and SySts requires both being able to correctly identify the 
origin of the ionised material and fully characterising the mor¬ 
phology and kinematics of the nebula. Here, two well known 
extreme bPNe (M 2-9 and Mz 3) and one confirmed symbiotic 
nebula (Hen 2-104) are investigated. It has been proposed that 
SySt mechanism has led to the formation of M 2-9 (|Doyle 


et al. 2000 [Livio & Soker 2001 [ Schmeja & Kimeswenger|2001 [ 


Smith & Gehrz |2005[ Corradi et al. |2011|l, Mz 3 ([Schnieja & 


Kimeswenger|2001[|Guerrero et al.|2604HSmith & Ciehrz|2005 , 


and H en 2-104 ([Whitelock 1987 Lutz et al. 1989 ISchwarz et ai.| 

|1989[[^rradi et al.|2001[|Santander-Garcfa et al. 120087! 

Previous studies on the morpho-kinematics of the structures 
of M 2-9 and Mz 3 have revealed similar morphological fea¬ 
tures (i.e., bipolar lobes), as well as the presence of material 
expanding from their central core up to a few hundred km s * 


([Lopez & Meabum|1983 Balick T989t|Solf 2000 Guerrero et al. 


2002| l. These materials have been proposed to be associated with 
highly expanding stellar wind from a cold companion. |Lykou| 


jet al.j ( |201 l| l and jCastro-Carrizo et ^ \20\2) found further ev¬ 
idence of a dusty circumstellar disk in the central region of M 
2-9. This result is consistent with the presence of a close binary 
system but has yet to be confirmed. Regarding Hen 2-104, its 


symbiotic nature has been confirmed by Whitelock & Munari 
( 1992|l after detecting a Mira companion in its central zone. The 


equatorial expansion velocity of this nebu la, which has a value of 
12 km s ' (Santander-Garcfa et al. 2008 1 , is significantly lower 
than those found in M 2-9 and Mz 3. 

We use optical and infrared diagnostic diagrams such as a 
T4363/H7 vs. d5007/Hj8 and a 2MASS NIR colour-colour di¬ 
agram to demonstrate that these methods discriminate well be¬ 
tween single ionising stars and symbiotic central stars. This tech¬ 
nique ( jGutierrez-Moreno et al. |1995| l demonstrates that large 
numbers of PNe and symbiotics (D- and S-types) can be read¬ 
ily separated. Thus, using these diagnostic diagrams, we charac¬ 
terise the central sources of M 2-9, Mz 3, and Hen 2-1 04. We 
then use the 3-D morpho-kinematic code called shap^ (Steffen 


et al. 201 l|l to carefully reproduce the structural and spectral fea- 


Garcia et al.||2007| have also been used to separate symbiotic 

nebulae from genuine PNe. On a J-H/H-Ks colour-colour dia¬ 
gram, the location of symbiotic nebulae can be found quite far 
from those of PNe and lie in the region of increasing J-H and 
H-Ks, as expected for D-type symbiotic nebulae (see Fig. 1 from 
|Phillips|2007| l. Nevertheless, the locus of symbiotic nebulae are 
also mixed with other dusty celestial objects such as Be, B[e], 
WR stars, or H ii regions, that make their classification unreli¬ 
able. Recently, an additional diagnostic diagram was proposed 
as a way to distingush symbiotic nebulae from genuine PNe by 
Corradi and collaborators, based on their r'-Ha and r'-i' colours 
( Corradi et al.|2008 2010[ Rodriguez-Flores et al.|2014) l. 

Interacting binaries have been shown to produce complex en¬ 
vironments, which makes them very useful for studying the late 
stages of stellar evolution, as well as the mass loss exchange 
between the components. It is generally accepted that accretion 
disks are needed to launch jets and possibly ansae (jGongalvesj 


tures of each nebula in three dimensions. For two of the three 
nebulae, it appears highly likely that a symbiotic central system 
is present and in all three, a complex multicomponent velocity 
structure is present within the overall bipolar morphologies. 

This paper is organised as follows: archival and source data 
and their use is presented in Sect.j^ Analyses and results of each 
topic are presented in Sect. This section includes: the mod¬ 
elling of optical kinematic data for the objects M 2-9, Mz 3, 
and Hen 2-104, using shape, and a spectroscopic analysis of 
M 2-9, Mz 3, Hen 2-104, and a number of other symbiotic 
and PNe sources in the optical and infrared using a /14363/H7 
vs. /15007/Hj6 diagnostic diagram and a 2MASS J-H/H-Ks dia¬ 
gram. Discussion and conclusions are found in Sects. |^and|^ 
respectively. 

2. Archival data 

2.1. Long-slit echelle spectra 

The observed long-slit echelle spectra used to assist with the 
creation of our model of M 2-9 were obtained from the 
San Pedr o Martir kinematic catalogue of Galactic Planetary 
Nebula ( [Lopez et al.|2012 i. Their observations were taken on 
June 2004 at the 2.1 m telescope in San Pedro Martir National 
Observatory, Mexico, with the Manchester Echelle Spectrometer 
(MES; [Meaburn et al. [[20031 1. Each slit had a slit width of 0.9" 
and a velocity resolution of 5.75 km s '. The spatial scale and 
sampling used for each slit were 0.312" pixel ' and 2.3 km s ' 
pixel ', respectively. The spectra in [N ii] were chosen for the 
morpho-kinematic modelling of M 2-9 because of the lower 
thermal broadening compared to the Ha line (thermal width of 


Ha and [N ii] at 10,000 K are 21.4 km s“' and 5.8 km s“', re¬ 
spectively). 

The observed long-slit echelle spectra used to assist with the 
creation of our model of Mz 3 were obtained from IGuerrero et akl 
( 2004| l. Their observations were taken on June 2002 using the 
echelle spectrograph on the CTIO 4 m telescope at Cerro Tololo, 
Chile. A slit width of 0.9" and a velocity resolution of 8 km s ' 
were used for each slit position. The spatial scale and sampling 
used were 0.26" pixeL' and 3.7 km s“' pixeL', respectively. 
Optical seeing was better than 1.2" during the observed period. 

The observed long-slit echelle spectra used to assist with the 
creation of our model of Hen 2-104 were obtained from ICorradil 


‘ SHAPE is accessible online at http://www.astrosen.unam.mx/shape/ 


http://kincatpn.astrosen.unam.mx 
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|et al.| ( |2001[ ). Their observations were taken on June 2000 with 
the CTIO 4 m telescope. The slit width used was 1.5" and the 
spatial scale along the slit was 0.54" pixel The sampling and 
resolution used were 3.7 km s“^ pixeH* and 13 km s“\ respec¬ 
tively. The seeing varied between 1.5" and 1.9" during the ob¬ 
servations. 


2.2. Archival data for diagnostic diagrams 


The emission lines used for the objects in the optical diagnos¬ 
tic diagram and hence the separation between SySts and PNe 
were obtained from various works in the literature (see Tables]^ 
and 1^ in the appendix). The 2MASS colour magnitudes of the 
same sample of SySts and PNe were obtained from the All-Sky 
Catalog of Point Sources NAS A/IPAC Infrared Science Archive: 
Cutri et ah] ( 2003a|b| l. All colour magnitude values are sum¬ 
marised in Tables |^and|2] 


3. Analyses and results 

3.1. Morpho-kinematic modelling of M 2-9 


The Galactic PN M 2-9, also called Minkowski’s Butterfly 
Nebula or PN GOlO.S-tlS.O ( [Acker et al.|1992 1 , is an enigmatic 
object with intense aspherical flows and displays a highly col¬ 
limated bipolar structure (120" x 12", [Schwarz et al.|1997| l. Its 
bipolar structure contains two separate shells, an inner bulb and 
an outer neutral shell ( [Corradi et al.|20lT] ), which are thought to 
be formed from two equatorial rings/tori ([Castro-Carrizo et~ar 


\2Q12) found at the nebula’s centre. The nebula also has a feath¬ 
ered component which extends all the way out to a pair of dusty 
blobs. Sub-structures in the form of arcs and surface brightened 
spots are also apparent. 

Apart from the dusty blobs and extended material (matter in 
the region beyond which the feathers disappear in the [N it] im¬ 
age), we have managed to construct and constrain the morpho- 
kinematics of M 2-9 and its components using a 3-D astrophys- 
ical modelling tool called shape (version 5.0). One aim involved 
creating the 3-D morphology resembling the nebula and its com¬ 
ponents using observed spectra (optical long-slit echelle spec¬ 
tra), which are compared and matched up with synthetic spectra 
generated from our model. It should be noted that the same pro¬ 
cedure described is also followed for the constmction of Mz 3 
and Hen 2-104. 

Development of the model involved the use of five individ¬ 
ual components; an inner closed-end bipolar component (bulb), 
a surrounding feathered shell (feathers), an outer bipolar com¬ 
ponent (neutral shell), and two tori (inner toms ring 1 and outer 
torus ring 2). The arcs (N2,S2) and surface brightened regions 
(N1,N3,S1,S3) where created by brightness enhancing of the 
bulb component surface. The names of the components were 
ado pted from[Kohoutek & Surdej[ ( [1980| ), [Corradi et al.[ ( [201 l[ l, 
and Castro-Carrizo et al. { 2Q12\ . Each individual component 
was modified using modifiers such as; density, velocity, shell, 
squeeze, rotation, etc. For example, the density modifier sets the 
spatial distribution of density, the velocity modifier applies a ve¬ 
locity field as a function of position, the shell modifier creates a 
shell of given thickness, the squeeze modifier compresses or ex¬ 
pands an object as a function of position along the local symme¬ 
try axis of the modifier, and the rotation modifier allows a chosen 
object (sphere, toms, etc,.) to be rotated about a specific point. 
A number of density and squeeze modifiers were customised in¬ 
teractively instead of analytically to carefully produce accurate 
results. 


Table 1. SHAPE components and velocities of M 2-9. 


Component 

location 

velocity (km s ') 

Bulb 

northern pole 

96 + 3 (27 + 3) 


southern pole 

40 + 3 (-29 + 3) 

Feathers 

northmax 

85 + 3(16 + 3) 


SOUthmax 

39 + 3 (-30 + 3) 

Neutral Shell 

northmax 

97 + 2 (28 + 2) 


SOUthmax 

41+2 (-28 + 2) 

'Ring 1 

— 

72.9 + 0.1 (3.9+ 0.1) 

' Ring 2 

— 

76.8 + 0.1 (7.8+ 0.1) 


Notes. The velocities shown are values with respect to Vsys = +69 km 
s“'. The values in parentheses are given as V- Kys (velocities corrected 
for Vsys). Here north^ax and south,^ax refers to the material found fur¬ 
thest out since these typ e of objects do not have a peak polar region, 
"velocities adopted from Castro-Carrizo et al. (2012 1 . 


In order to mimic the observed spectra and constrain the 
kinematics of M 2-9, each stmcture and sub-stmcture required 
a velocity component. The velocity function in shape is given by 
I f I (r) = A H- B(r/ro) + C(r/rf)^, where r is the initial arbitrary 
(constant) radius given to the structure, whereas all the other pa¬ 
rameters are variables. The pq parameter is used to give the rate 
of change of velocity with respect to the radius r. Homologous 
and radial vector fields were also assigned to the different com¬ 
ponents to account for their type of expansion. There is some 
evidence however of localised deviations from a radial flow in 
the outer lobes of M 2-9, where the addition of a modest lat¬ 
erally expanding velocity component provides a better fit. See 
Table. for the velocities of each structure and the appendix for 
the 3-D mesh images of M 2-9. 

The four long-slit spectra of [Lopez et al. (2012 1 were used 
to create and test our model of M 2-9. All four slits are posi¬ 
tioned parallel to the nebula’s symmetry axis (north to south). 
Slit c is placed directly down the symmetry axis, whereas slit a, 
b, and d are offset from the axis by -5.0", -3.0", and -1-3.5", re¬ 
spectively. It should be noted that the observed long-slit spectra 
from [Solf[ ( [2000[ l and [Smith & Geh^ ( [2005| l were also used to 
help assist with the modelling. A slit width and slit length of 1" 
and 60" respectively were used to render the synthetic spectra. 
In each case, the synthetic spectra are rendered at the same ve¬ 
locity range (330 km s“*) and spectral resolution (5.75 km s“') 
as the observed spectra. We should also note that the velocities 
used for the modelling of ring 1 and ring 2 were adopted from 
Castro-Carrizo et al. (2012 1 due to the intensity of [N n] emission 
from the central region (see observed spectra b and c in Fig.[2l. 
The systemic veloci ty Vsys for the ob served and synthetic spectra 
is -1-69 km s“' (see Smith & Gehrz|[2005 1. Both sets of spectra 
are presented alongside each other in Fig. ^ The HST difference 
imag^and rendered shape image can be seen in Fig. 

Our best-fit model (including all bipolar components) has an 
inclination of its symmetry axis with respect to the plane of the 
sky of i - 17°, which agrees with the value found by [Zweigl^ 
et al. (19971 and Solf[ ( 2000|, but is slig htly smaller than the 
value (i - 18°) found by Corradi et al. (20111. A marginally 
smaller value of i - 


16° was found by 
whereas values of i - 15° ([Schwarz et a. 


([Goodrich|1991[l were also determined. 


Lykou et al.[ ([201IJ, 
.[[1997| l and i = 11° 


http://apod.nasa.gov/apod/ap 130915.html 
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Vhkl [km s’*] 


Fig. 1. Observed [N ii] /16583 (top row) and synthetic (bottom row) spectra for M 2-9. The brightness enhanced features are labelled 
on both sets of spectra for clarity. The systemic velocity Vsys = 69 km s“'. The observed spectra shown are from Lopez et al. (2012 i. 



II 


Fig. 2. Grayscale NST [N ii] image (left) and shape model (right) of M 2-9. The main components and features of the nebula are 
labelled in the [N ii] image. Both equatorial rings seen across the minor axis in the shape image are best observed in the '^CO J - 2-1 
line emission ([Castro-Carrizo et al.|2012[). The height and width of each box is 48"x 48". 


The distance to M 2-9 has been routinely quoted throughout 


the literature and is said to have values of 650 pc (Schwarz et al. 

|19971|Doyle et al.|2000 |Solf|2000l |Castro-Carrizo et al. 

20121, 

1.2 kpc (|Lykou et al. 201 l|l, and 1.3 kpc (|Corradi et al. 

20111. 


In this study of M 2-9 we adopt the distance of 650 pc. At this 
distance we find a kinematical age of 915 + 90 yr for the bulb. 


Castro-Carrizo et ^ ( 2012| l found a similar age for the inner 
molecular torus (ring 1) of ~900 yr at 650 pc, which suggests 
that the bulb may have formed just after the formation of ring 1. 
The kinematical age we find for the neutral shell is 1430 + 90 
yr, whereas Solf ( 2000| l determined an age of 1300 yr. Castro- 
ICarrizo et al.|P012| l and [Zweigle et al.| ( |1997| l found values of 
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1400 yr and 1365 yr for ring 2, respectively, which suggests, like 
that of the bulb, that the neutral shell could have formed after the 
formation of ring 2. This evidently indicates that both rings are 
the equatorial counterparts of these bipolar outflows. 

Regarding the feathered component, we were unable to de¬ 
termine an accurate kinematical age because our modelled spec¬ 
tra does not include the extended material and dusty blobs (N4 
and S4) which most likely makes part of the feathers. |Corradi| 
et al. (20111 found the kinematical age of these blobs to be 
2500 yr for a distance of 1.3 kpc and a velocity of 147 km s '. 
However, considering a distance of 650 pc with the same veloc¬ 
ity of 147 km s“*, the kinematical age of the blobs would be 
~1250 yr. This indicates that the neutral shell was the earliest 
ejection event, then followed by the feathers and dusty blobs, 
and most recently the bulb. 

3.2. Morpho-kinematic modeiiing of Mz 3 

Menzel 3 (PN G331.7-01.0), or Mz 3 as it is commonly known, 
is an intrinsically bright young bipolar PN, which has arguably 
the most complex bipolar morphology seen to date. Mz 3 con¬ 
sists of three nested pairs of bipolar lobes {BLl, BL2, and BL3), 
an equatorial ellipse-like feature {EE2), and three knots (knot 1, 
knot 2, and knot 3), which we adopted from ([Guerrero et al. 


Table 2. SHAPE components and velocities of Mz 3. 


2004 hereafter Gu04). BLl appears to have a closed-lobe mor¬ 


phology, whereas BL2 and BL3 have open lobes with cylindrical 
and conical shapes, respectively. The faint equatorial ellipse of 
emission appears to be aligned along the equator of the nested 
bipolar lobes (see Fig. 1 of Gu04 for a more noticeable view of 
this component). The three faint knots, two (knot 1 and knot 2) 
of which are located just above the apex of the northern lobe of 
BLl, and one (knot 3) found just below the southern lobe, were 
initially discovered by Gu04. 

To get a better understanding of Mz 3’s complex morphol¬ 
ogy and kinematics, a shape model of the nebula was created. 
Completion of this model included seven components, two of 
which contained a number of sub-components. For example, the 
inner bilobed structure BLl, or IBL in jRedman et al.jpOOOj ), con¬ 
sists of two inner bubbles (Bubble 1 and Bubble 2) that represent 
what appears to be a pair of inner bilobes, and the top region 
of each lobe of BLl contains twelve sub-bubbles that represent 
the foamy texture, or blister-like structures in Gu04, whereas the 
outer bipolar conical structure BL3, or rays in ( jSantander-Garcfa 
et al.|20()4 hereafter SG04), is made up of fifty-eight individual 
ray-like structures. See Table|^for all components and velocities 
used to create our model of Mz 3. 

The observed spectra (Fig.|^ used to assist in the creation of 
Mz 3 are those from Gu04. Their high-dispersion spectroscopic 
observations allowed us to determine kinematic properties of the 
different components, as well as tightly constrain the nebula’s 
morphology. The four slit positions at P.A.s -28°, 8°, 52°, and 
98° were considered because of the important regions they cover 
on the nebula (see Fig. 3 from Gu04 for a detailed image of 
these slit positions). The different morphological components of 
Mz 3 are labelled on the observed spectra, as clearly shown in 
Fig. 0 The height and velocity range for both sets of spectra are 
100 arcsec and 600 km s ', respectively. The difference imag^ 
which is a composite grayscale HST image, and the rendered 
SHAPE image are shown in Fig.|^ See appendix for the 3-D mesh 
images of Mz 3. 

All three pairs of bipolar lobes share the same axis of sym¬ 
metry (Gu04, SG04) and each one is modelled with an inclina- 


Component 

location 

velocity (km s ') 

BLl 

northern pole 

100 + 5 


southern pole 

-80 + 5 

BL2 

northniax 

110 + 7 


SOUthmax 

-100 + 7 

BL3 

northmax 

140 + 10 


SOUthmax 

-100+10 

EE 

equatorial 

200 + 10 

Knot 1 

— 

30 + 2 

Knot 2 

— 

40 + 2 

Knot 3 

— 

-30 + 2 


Notes. The velocities shown are given as V - i4ys- 


tion of its symmetry axis with respect to the plane of the sky of 
i — 12°. An inclination between 12° and 20° was found by Gu04, 
whereas SG04 find a value between 14° and 19°, and Lykou et al. 
( 2007| l proposed a value of i = 16°. 


The distance to Mz 3 remains uncertain and a range of values 


have been previously published; 1 kpc (Lopez & Meaburn 

1983 1, 

1.3 kpc (jCahn et al. 

1992 

iPottasch & Surendrranath 

20051, 

1.4 kpc ([Lykou et al. 

2007 

1, 1.8 kpc (Cohen et al.[[1978i, and 

2.7 kpc ([Kingsburgh & English|1992 1 . Due to this uncertainty, 


http://www.astro.washington.edu/users/balick/WFPC2/index.html 


we will show explicitly the dependence of the kinematical age 
on distance, like those shown in Gu04 and SG04. 

SG04 show BLl to have a kinematical age of fD ' (yr kpc *) 
= 670 yr with a possible range of 550-710 yr, whereas Gu04 find 
BLl to be 600+50(Z)/kpc) yr, where D is the distance in kilopar- 
secs. We find a kinematical age of 650 + 35(D/kpc) yr. Also, the 
inner bipolar lobes (bubble 1 and bubble 2) found within BLl 
interestingly share the same kinematical age as BLl for a Vmax = 
65 km s ' and ~9" displacement from the nebular centre. This 
may suggest that the inner bipolar lobes and BLl most likely 
formed at the same time. 

The knots, found just above the poles of BLl, were modelled 
with displacements (values are from Gu04) of 19" (knot 1), 24" 
(knot 2), and 19" (knot 3) from the central star and along the 
symmetry axis. The knots are shown to have kinematical ages of 
1020 + 70(Z)/kpc) yr (knot 1), 1200 + 60(D/kpc) yr (knot 2), and 
1020 + 70(D/kpc) yr (knot 3), taking into consideration that they 
lie along the symmetry axis of BLl and follow that of a Hubble 
flow. 

The morphology of BL2, or columns in SG04, is quite dif¬ 
ferent to that of BLL It has a highly collimated cylindrical shell, 
which increases gradually in diameter (d) along each lobe from 
the end point of the curvature (d ~2Q") near its centre to its 
furthest observable point out (d ~30"). BL2 also has a tight 
waist (2.8+0.2") and a relatively thick shell (2.5+0.2"). We find 
a kinematical age for BL2 to be 1000 + 50(Z)/kpc) yr, which 
confirms the same age found by Gu04. A value of 875 yr, with a 
possible range 830-1000 yr , was acquired by SG04 for BL2. 

A more complicated structure to create was that of BL3. 
Although it is a single structure by nature, here it is made up 
of fifty-eight individual rays to give its stringy-like appearance. 
Each ray was created as a long thin cone to give the effect of 
expansion of the material. These rays are positioned in such a 
way that they give a conical bipolar shape to the entire structure. 
The opening angles for the lobes of BL3 are ~55° (47° and 50° 
in SG04 and Gu04 respectively) with respect to the axis of sym- 
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(V-Vsys) [km s-'l 

Fig. 3. Observed [N ii] /16583 (top row) and synthetic (bottom row) spectra for Mz 3. The RA.s and nebular components are labelled 
on the observed spectra. The velocities are corrected for the systemic velocity of the nebula. The observed spectra shown are from 
[Guerrero et al.| ( |2()04| l. 



Fig. 4. Grayscale composite HST image (left) and shape model (right) of Mz 3. The main components of the nebula are shown in 
the HST image. The equatorial ellipse (EE) is difficult to observe in the composite HST image but is best seen in Eig. 1 of Gu04. 
The height and width of each box is 96"x 96". 
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metry, and are positioned around a 360° angle with the nebular 
centre as the point of emanation. For a maximum extension of 
the rays ~55" (60" in Gu04) and adopting a Hubble-like flow, 
we determine a kinematical age of 1800 + 100(D/kpc) yr, which 
agrees with the value found in Gu04. A value of 1600 yr, with a 
possible range 1400-1800 yr, was obtained by SG04 for BL3. 

On the other hand, the EE component was modelled as an 
elliptical shell with a higher expansion velocity at its equatorial 
regions than at its poles. Like that of the three bipolar compo¬ 
nents, we modelled EE with an inclination of 12° but with its 
symmetry axis (that connecting its poles) tilted by -i-8° (RA. 0°). 
Using our modelled parameters, we hnd a kinematical age of 
950 + 50(D/kpc) yr, which is in close agreement with the value 
found by Gu04 (1000 + 50(Z)/kpc) yr). Therefore, considering 
the kinematical ages derived, and those compared with Gu04 
and SG04, the rays of BL3 would be hrst produced, then fol¬ 
lowed by the combination of BL2, EE and the knots, and most 
recently the BLl structure. 


3.3. Morpho-kinematic modeiiing of Hen 2-104 

The Southern Crab (PN G315.4H-09.4), or Hen 2-104, is a sym¬ 
biotic nebula that exhibits an ionised bipolar outflow as a re¬ 
sult of a symbiotic binary system at its centre. Hen 2-104 was 
once classihed as a PN until Whitelock] ( 1987| l recognised a 
long-period Mira at its centre. The formation of bipolar nebu¬ 
lae around symbiotic stars are said to be the result of complex 
interactions, and Hen 2-104 clearly appears to be one such ex¬ 
ample. 

The nebula consists of three main components: an inner 
hourglass, an outer open-ended hourglass, and a pair of po¬ 
lar jets found along the symmetry axis of the system (see 


Corradi & Schwarz) 1993 [[Corradi et al.|2001[|Santander-Garcfa] 


et al. 20081. The outer bipolar nebula has a span of 70"x 35 
(Santander-Garcia et al.|[2008l in the plane of the sky, and the 


inner hourglass is about five times smaller in scale (Corradi 
et al.||200T l. The inner and outer hourglass structures share the 
same symmetry axis (Corradi et al.|200l] Santander-Garcia et al. 


2008|l and are inclined by 55° respect to the line of sight, with 


the southern lobe pointed towards us. Corradi et al. (|2001| l found 
an inclination i - 60°, whereas [Santander-Garcia et al.| ( [2008| 
determined a value i = 58°. Three symmetrical rings are seen 
on the higher latitudes of each lobe of the inner hourglass, as 
shown by [Corradi et al. pOOlj l. Where the inner hourglass ap¬ 
pears intact, the outer hourglass is highly fragmented in the form 
of knots. 

The overall structure of Hen 2-104 has been modelled be¬ 
fore by Corradi et al. (20011 and Santander-Garcia et al. ( 2008) 1, 
both of whom performed a comprehensive study of the nebula’s 
inner and outer hourglass’ and its high-velocity polar knotty out¬ 
flow. However, in order to hnd any potential systematics in the 
distribution of the small-scale structure, we are interested in the 
full 3-D distribution of the knots and hlaments. Eor our morpho- 
kinematic analysis of Hen 2-104 we used the observed long-slit 
spectra (Pig. [^ from Corradi et al. (2001 1 . Six long-slit posi¬ 
tions with PA.s from h- 122° to -t-272°, in increments of -1-30° and 
through the nebular centre, were used to help create our model 
of Hen 2-104. The slits positioned at PA.s -1-122° and -1-212° 
pass directly through the major axis (axis of symmetry) and mi¬ 
nor axis, respectively. The HST difference imag^and rendered 
SHAPE image are shown in Pig. [^ See Table [^ for a list of the 
components used to create our model of Hen 2-104. 


http://www.spacetelescope.org/images/opo9932a/ 


Tables. SHAPE components, edge distances, and constants of pro¬ 
portionality for the model of Hen 2-104. 


Component 

edgGinner ( ) 

^^S^outer ( ) 

k (km s V”) 

Inner bipolar 

— 

5.6 

6 

Knots (red) 

10.3 

23.5 

2.5 

Knots (blue) 

12.3 

20.0 

2.5 

'Bipolarsmooth 

— 

— 

2.5 

'^Jet knoti-ed 

— 

— 

10 

'■'■Jet knotbiue 

— 

— 

— 


Notes. The velocity law is v(r) = k(r/ro), where k (given in 
km s“* arcsec^') is the constant of proportionality that maps the po¬ 
sition to velocity, ro is a fiducial reference distance, and r is the distance 
from the centre as measured in arcsecs. The deprojected distances to 
the edges of the components are given in arcsecs. ' The outer edge of 
the smooth bipolar component ends where the knots begin. " The Jet 
knotred has a velocity of ~130 km s“* at 13", whereas the velocity of Jet 
knotdiue (located at 11") is unconstrained, since it is not covered by the 
spectroscopy. 


The highly knotty and filamentary structure of Hen 2-104 re¬ 
quired a slightly different modelling approach compared to the 
other objects. We therefore first produced a cylindrically sym¬ 
metric mesh model of the overall structure, taking into account 
the HST image and observed long-slit spectra. In particular, we 
considered the inner smooth low level emission and the knots 
and filaments further out. It turned out that the clear transition 
from a smooth emission distribution to the knots and filaments at 
an angular distance from the centre of about 17" was a very use¬ 
ful constraint for the model. In the regions of ambiguous posi¬ 
tioning in the image (slit positions h- 182° and h- 242°), where the 
line of sight is nearly tangential to the lobes, we restricted the po¬ 
sitioning of the knots to a region (outside the red and blue rings 
shown in the mesh images of Hen 2-104 in the appendix) where 
they start in the locations that are not ambiguous. Purthermore, 
this constraint helped to identify two inner knots near the pro¬ 
jected axis that we now attribute to the collimated outflow along 
the axis, which are indicated in the bottom two mesh images as 
‘Jet’. 

To model the detailed position and structure of the individ¬ 
ual knots we applied particles to the positions of the knots and 
filaments, shape incorporates a special 3-D cursor tool whose po¬ 
sition can be limited to a surface mesh and then located to any 
position with the help of the mouse or a pen device on a graphics 
tablet. The rendering size and brightness of the particles can also 
be controlled. Bright knots have been given a higher brightness 
than lower surface brightness filaments. This brightness assig¬ 
nation is not meant to reproduce the actual brightness in detail 
but is meant to be a rough indication of the relative brightness 
of the features, which is sufficient for the current purpose of our 
analysis. 

The accuracy of the positioning of the particles in the re¬ 
gions where the lobe surface is seen nearly face on is close to 
the level of the HST image resolution. In regions where the line 
of sight is nearly tangential to the lobe surface, the positioning 
of the knots and filaments relies mainly on the position in the 
lower-resolution long-slit spectra and the additional constraints 
mentioned above. It turns out that the positioning is still very 
accurate, probably only a factor of two less accurate than in the 
other regions. 
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(V-Vsys) (km s-') 

Fig. 5. Observed [N ii] CTIO (top row) and synthetic (bottom row) spectra for Hen 2-104. Shown are the spectra at P.A.s. 122°, 
182°, and 242°, for simplicity. The radial velocities are corrected for the systemic velocity of the nebula. The observed spectra are 
from [Corradi et al.|P001| l . 



Fig. 6. HST [N ii] image (left) and shape model (right) of Hen 2-104. The shape model shows the small-scale structures (fragments) 
of the outer shell in Doppler colour to display their direction of motion. The height and width of each box is 62.5"x 62.5". 


3.4. The optical PNe versus SySts diagram 

In the last twenty years, several efforts have been made to sepa¬ 
rate symbiotic-type nebulae from genuine PNe, mainly based on 


observations in the optical regime. One effort of this kind was 
the BPT diagram ( Baldwin et al.|198T| l that attempts to separate 
extragalactic emission-line objects according to their excitation 
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mechanism. Here we present a diagnostic diagram based on the 
/14363/Hy vs. /15007/Hy8 emission line ratios for the separation of 
genuine PNe from SySts, as proposed by|Gutierrez-Moreno et al. 
( |T995l l. Fig. |7] plots d4363/Hy vs. d5007/H/3 for the three ob¬ 
jects studied here (M 2-9, Mz 3, and Hen 2-104) together with 
a signihcant sample of known SySts and PNe. Indeed, SySts and 
PNe are well separated in two different loci. Only two objects 
are found to be misplaced in the region of young PNe (Hen 2- 
171 and H 1-36). In particular, H 1-36 is a well known D- 
type SySt with a variable Mira star (Medina Tanco & Steiner 


1995|l. The misplacement of known SySts in the optical diagnos¬ 


tic diagram has previously been reported for other known SySts, 
such as Hen 2-147 and the extragalactic SySt (ICIO SySt-1); 
in addition of two new likely SySts found recently in the dwarf 
spheroidal galaxy NGC 205 ( Gongalves et al.|2014| l. [Gutierrez - 
Moreno et al. (|1995 i have also mentioned that the separation of 


D-type SySts and young PNe is not always possible with this 
diagnostic diagram because of the high density nature of these 
two types of sources. The bifurcation of the data points in this 
diagram is the outcome of the different densities between the 
two groups of objects. For instance, SySts have a density several 
times higher than PNe and this results in higher attenuation of 
[O III] 5007 A emission to that of [O m] 4363 A emission. This 
is the reason for the emission and hence the clear separation of 
the SySts from the PNe in the optical diagnostic diagram. 

Regarding M 2-9 and Hen 2-104, these objects are well lo¬ 
cated in the regime of SySts, providing additional evidence for 
their symbiotic nature, whereas Mz 3 is placed in the bottom- 
left corner of the diagram, making its classihcation more uncer¬ 
tain. In particular, the emission-line ratios of Mz 3 have been 
estimated, se parately, for i ts inner region and also for the outer 
bipolar lobes (Smith 2003 i. The A4363 A line is found stronger 
in the inner region compared to the bipolar lobes, resulting in a 
higher /14363/Hy, which moves Mz 3 deeper into the regime of 
SySts (see Fig. |^. However, its real nature still remains debat¬ 
able. On the other hand. Hen 2-104 is a known D-type SySt with 
a Mira variable companion ( Whitelock & Munari|1992 1 . As for 
M 2-9, several authors have classihed this enigmatic object as 
a SySt but its true nature still remains debatable. Nevertheless, 
the presence of two compact circumstellar disks in M 2-9 have 
recently been conhrmed ( Lykou et al.|201 1[ Castro-Carrizo et al. 
|2012^ , providing further evidence of a likely binary system in its 
core, which is consistent with the SySt scenario, as much as with 
a genuine PNe with a binary central star. 

Another PN once thought to be a SySt is MyCn 18 due to 
the similar characteristics it has with symbiotics. O’Connor et al. 
( 2000[ l; [Soker & Kashi ( 2012| l and Clyne et al. ( |2014| l suggested 
the possibility that MyCn 18 might have been a SySt, however, it 
is clearly placed in the PNe regime of the optical diagnostic di¬ 
agram (Fig.[^. This further suggests that MyCn 18 is a genuine 
PN. 

Also of importance are the uncertainties of the emission lines 
([O III] 4363 A, [O III] 5007 A, H/3, and Hy), which are lower 
than 15% for most cases. However, for some PNe the uncertainty 
of the weak [O iii] 4363 A line can be between 20^0%. This 
makes objects located in the small region at the bottom left-hand 
corner (0 < /14363/Hy < 2.5 and 0 < A5Q01/H/3 < 1), and those 
bordering between regions, uncertain for proper identihcation. 
Mz 3 is located in the small boxed region, which makes identi¬ 
fying its symbiotic nature problematic. One additional cause that 
may affect the value of the 4363/Hy line ratio, and thus the clas¬ 
sihcation of the sources, is the Hg 4359 A tulleric line, which is 
sometimes easily mistaken for the faint [O iii] 4363 A line. 


3.5. The 2MASS PNe versus SySts diagram 

Unlike the optical diagram just discussed, in the 2MASS J-H/H- 
Ks diagnostic diagram the regime of SySts is mixed with other 
dusty celestial objects like Be, B[e], WR stars, as well as young 
dusty PNe and H ii regions, making the identihcation of these 
objects difficult. Fig. also represents the 2MASS J-H/H-Ks 
diagram for the same sample of SySts and PNe as in the optical 
plot just above it. This 2MASS diagram clearly shows that SySts 
dehne a specihc region where all of them are placed. Only one 
SySt is found to be misplaced, namely SS 122. This well known 
SySt is, however, found to be well placed in the regime of SySts 
in the optical diagram. We should also mention here that two 
SySts (Hen 2-171, H 1-36), which are found to lie in the region 
of young-PNe in the optical diagram, are also well placed in the 
regime of SySts in this infrared diagram. 

According to their near-infrared colours and spectral energy 
distributions (SEDs), SySts are also classihed into three sub¬ 
groups; S- (Stellar), D- (Dusty), and D'-type. The S-type SySts 
show NIR temperatures in the range 2500-3000 K, which is con¬ 
sistent with the effective temperature of a K- or M-type giant, 
whereas D-type temperatures of 1000 K indicates a compan¬ 
ion obscured within a warm dust shell (Mikolajewska 1997| l. 
For comparison, we also plot the NIR colour values for a group 
of K and M giants ( van Belle et al.||T999| l. These objects share 
the same regime with the S-type SySts, which is expected with 
a smaller range of J-H values. The difference between D- and 
D'-type SySts is the peak in their SED’s, with the former and 
latter peaking at 5-15 pm and 15-30 pm respectively pvison 
et al.|1995) . On the other hand, the group of PNe covers a larger 
range of H-Ks colours from approximately -0.7 to 2.0, whilst 
the majority are concenPated in an even smaller region with J- 
H colours between -0.5 and 0.5. There are, however, few PNe 
lying in the region of SySts. 

M 2-9, Mz 3, and Hen 2-104, are shown as red pentagons in 
the infrared diagram, and lie within the region of D-type SySts. 
This is consistent with the Mira companion found in Hen 2-104 
( [Whitelock & Munari|| 19921 , but, thus far, there is no observa¬ 
tional evidence of a Mira companion in either M 2-9 or Mz 3. 
The SEDs of these two objects indicates the presence of two 
stellar components, a hot and cold star (Smith 2003 [ Lykou et al.| 
20111. In particular, the recently modelled SED of M 2-9, con- 
sPucted by Lykou et al. ( |201 1 1 using various photometric and 
spectroscopic data, shows a peak at ~19 pm, which is similar 
to those of D-type SySts. Nevertheless, this could simply prove 
the presence of a dense circumstellar envelope, not necessarily 
indicating that it is a SySt. Regarding the PN MyCn 18, it is 
found to be placed within the region where PNe are located in 
the 2MASS diagram, and if taken in conjunction with the optical 
diagram, only further implies its PN nature. 


3.6. Ha broadening with doubie-peak prof ties 

One particular velocity component that we did not reproduce in 
our models is the high, broadened velocity that emanates from 
the central region of each object (see the observed specPa in 
Pigs. BUI and|^. However, we thought it would be of interest to 
analyse the velocity profiles that they exhibit. Interestingly, we 
found that the Ha specPa from slit b and c of M 2-9, show a 
bright double-peaked profile, with the blue peak fainter than the 
red ( [Lopez et al.||2012 l. This kind of Ha profile is commonly 
found in SySts, as are broad single peak profiles ( [Van Winckel 
et al. 1 1993 1 Ivison et al. 1994 1 , and can be interpreted as the 
result of a sPong (high-mass loss rate) and fast stellar wind ex- 
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Fig. 7. /14363/Hy vs. /15007/Hy6 (top) and 2MASS NIR colour-colour (bottom) diagnostic diagrams of PNe and SySts. The triangles 
and squares correspond to PNe and SySts, respectively. The red and green pentagons represent the bipolar nebulae, M 2-9, Mz 3, 
and Hen 2-104, with MyCn 18 as a special case. The colours represent the different sources (see Tables.]^ andfrom which the 
data were obtained. The dashed box in the bottom left-hand corner of the optical diagnostic diagram defines a small region where 
the identification of the sources as SySts, Young PNe, or PNe is difficult due to the uncertainty in the weak [O iii] 4363 A line, and 
the two diagonal lines correspond to an [O iii] 5007/4364 line ratio of 13.1 (upper line) and 27.4 (lower line). The red and green 
crosses in the 2MASS NIR diagram correspond to K and M giants respectively, and the black arrow corresponds to 3 mag extinction 
in the V band. 
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Fig. 8. Ha line profiles of two known SySts, AS 245 and BX 
Mon, and two candidate SySts, M 2-9 and M 1-91. The data 
were obtained from the San Pedro Martir observatory (Akras 
et al. in prep). The x-axis represents the expansion velocities, 
where arbitrarily it has been given a zero velocity to the absorp¬ 
tion dips for direct comparison. 


pelled by the cold companion. The exact mechanism, however, 
remains undetermined. 

High-dispersion Ha line prohles of two known SySts, 
BX Mon and AS 245, and two candidate symbiotic nebulae, 
M 1-91 and M 2-9, are presented in Fig.[^ M 2-9 exhibits an ab¬ 
sorption feature similar to what we hnd for the two known SySts, 
whereas it appears less absorbed in M 1-91, indicating a differ¬ 
ent degree of obscuration along our line of sight (see |Torres-| 
Peimbert et al.|20f0 1 . Moreover, it may also be associated with 


the source type (S- or D-type; Van Winckel et al.|1993 1 or pos¬ 
sible Ha variations ( |Ivison et al. 11994) 1. The detection of the He 
II 6545 A Raman scattered line in M 2-9 by ~ 


Lee et al. 


(20011 in 


conjunction with the non-detection of any He recombination line 
suggests a highly obscured object. Moreover, M 2-9 and M 1- 
91 show signihcantly lower velocities (31 and 29 km s ') com- 
pai'ed to those of BX Mon and AS 245 (42 and 48 km s“'), which 
is probably associated with the different evolutionary stages of 
the binary system and/or the presence of an extended Raman- 
scattering origin ( |Lee et al.|2001| l. It should also be mentioned 
here that the observed [O iii] 5007/HjS and [O iii] 4363/Hy line 
ratios of M 1-91 (Rodrfguez et al. 200 1[ ) place its core in the 
regime of SySts in the optical diagnostic diagram, whereas the 
bipolar collimated lobes are placed in the region of young PNe, 
just like that of M 2-9. 

In addition to this analysis, and using the SPM catalogue of 
galactic planetary nebulae ( |L6pez et al. 2012| l, we revised all 
the available high-dispersion echelle spectra for the known PNe 
that posses a binary system at their core. Two out of 26 PNe 
were found to exhibit Ha double peak prohles with high ve¬ 
locities, namely Hen 2-171 and Sh 2-71. Surprisingly, Hen 2- 
171 is an already known SySt star, whereas Sh 2-71 is classi- 
hed as a PN. Sh 2-71 can not be applied to the optical diag¬ 
nostic diagram since there is no available emission line spectra 
from its core, however, its NIR colours place it hrmly within the 
regime of SySts in the 2MASS diagram. Recently, [Mocnik et al. 
(|2015|l conhrm the binarity of Sh 2-71’s central star, namely 


CSl, where the two companions are found to be a hot O-type 
subdwarf and a Be star. 

4. Discussion 

The long-slit echelle spectra provided us with the necessary in¬ 
formation about the 3-D structures of M 2-9, Mz 3, and Hen 2- 
104, allowing for the generation of self-consistent 3-D mod¬ 
els. The synthetic spectra generated from each model compares 
closely with their corresponding long-slit echelle spectra, which 
emphasises the accuracy of our models. 


4.1. The morpho-kinematics of M 2-9. 

What makes M 2-9 intriguing is not only its detailed morphol¬ 
ogy but the peculiar behaviour within its internal region. The 
motion of the features N1 and SI were initially thought to be 
knots moving in a lateral direction along the walls of the bulb, 
however, [Corradi et al. ( 2011| l described this apparent motion 
as excitation of the gas on the walls by a collimated spray of 
high velocity particles (jet) from the central source. In partic¬ 
ular, [Corradi et al.| ( |201 Ij ) described that the walls of the inner 
cavity (bulb) are ionised by shocks produced by the impact of 
the jet particles. This would explain for the lighthouse effect 
seen in the bulb. Livio & Soker ( 2001| l, on the other hand, pro¬ 
posed a beam of ionising photons (light beam). A light beam 
would require most of the illuminating source to be obscured, 
and would represent an additional source of localised ionisa¬ 
tion or scattered light. We accepted the former explanation by 
Corradi et al. (201 l|l as a more probable scenario for this un¬ 


usual motion, which is why our model of M 2-9 includes sur¬ 
face brightness enhancements for N1 and S1 instead of physical 
components (like knots). Closely observing the shape of these 
features on the HST [N ii] image, especially that of Nl, it can be 
seen that they follow the curvature of the wall, which strongly 
suggests that they are the result of illumination by excitation. 
Modelling of these features (and those of N2, S2, N3, and S3) 
worked well using this approach, as shown in the synthetic spec¬ 
tra of Fig. [T] 

It is also worth noting that the [N ii] emission from the cen¬ 
tral region is not reproduced in our model of M 2-9. This emis¬ 
sion is the result of the gaseous interactions between the two 
binary stars at its core, with the hot companion ionising the stel¬ 
lar wind from the cold companion. This fast expanding ionised 
wind is responsible for the broad wings in the Ha line, com¬ 
monly found in SySts. The emission of M 2-9 has been found 
by Trammell et al. ( |1995| l to be the combination of intrinsic/local 
emission and that of scattered emission from the central region. 
This scattered emission adds additional broadening to the emis¬ 
sion lines for the N2 feature. This implies that we may slightly 
overestimate the real expansion velocities of the components. 
However, our shape model of M 2-9 provides a very good inter¬ 
pretation of the current available data. 


On another note, an energetic rotating jet (Corradi et al. 
[201 1 1 is most likely due to an existing binary system, which 
could explain for the tori or rings seen in Castro-Carrizo et al. 
( 2012| l. This is why our model includes equatorial ring-like 
structures. Both equatorial rings appear to ht the waists of their 
corresponding coaxial lobes (see shape image of Fig. |^, which 
might explain for the formation of the bulb and neutral shell. 
[Castro-Carrizo et^ ( [2012| l also describes this system of rings 
as being ejected by one star, rather than a very late AGB or 
very early post-AGB phase, because of the small differences in 
their central offsets (0.34") and systemic velocities (0.6 km s“'). 
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Therefore, the nebula would have been formed in, at least, two 
separate events. Also, the observed [O iii] /14363/Hy and [O m] 
/15007/HyS emission line ratios from its central region indicate the 
presence of a symbiotic binary system (see Fig. |^. Moreover, 
the double-peaked Ha profile from its central region resembles 
those from other known SySts (Fig. [ 8 1 , which provides further 
evidence of this scenario (|Torres-Peimbert et al.|201()ll. See also 


Van Winckel et al. (1993 1 andflvison et al. (1994 1 for a descrip¬ 


tion of this. 

Furthermore, the formation of the feathers is not so easy to 
explain. The extended emission appears to originate along the 
walls of the bulb, which suggests the possibility that it is a sepa¬ 
rate event that formed sometime after the bulb itself. Moreover, 
if the bulb formed after the feathers then its outflowing material 
could have interacted and moved along the walls of the feathers 
( Kwok|1982[ ), which would explain their similar appearance. A 
pair of high-velocity dust blobs (N4 and S4) can also be found 
far from the poles of the bulb along the nebula’s axis of symme¬ 
try (see Fig. 1 from Corradi et al.|2011) and has extended emis¬ 
sion reaching all the way out to them, therefore, the feathers and 
dust blobs possibly formed at the same time. We should also note 
that the dust blobs and the extended emission connecting them 
were not included in our model due to the unavailability of long- 
slit spectral data for these components. Finally, as mentioned in 
Sect. im our SHAPE fit for the outer lobes supports a scenario 
in which an initially radial expansion becomes redirected, pre¬ 
sumably by dense equatorial material into a collimated, laterally 
expanding, outflow directed along the main symmetry axis of 
the nebula. This kind of recollimation is seen in outflows in star 
formation (Rawlings et al.|20T3 1. 


4.2. The morpho-kinematics of Mz 3. 

Modelling the components of Mz 3 proved a challenge due to 
their highly complex morphologies. The bilobed shell BLl, its 
internal frothy textures (bubbles and sub-bubbles), and the pro¬ 
trusions seen at the tips of the shell were all carefully con¬ 
structed, both morphology and kinematically, to match those 
seen in the observed high-quality long-slit spectra in Fig. and 
in the HST image of Fig. To replicate these observed spectra 
from Gu04 a second pair of inner bilobes (bubble 1 and bubble 
2) were included. This might suggest a relatively young ejec¬ 
tion event due to the presence of a binary system at its cen¬ 
tre. Interactions of this possible outflow with the apexes of BLl 
might explain for the interrupted blistered regions and protru¬ 
sions seen in the HST image of Fig. however, we find that 
BLl and these inner lobes share the kinematical age, which de¬ 
viates from this hypothesis. 

Also of interest are the knots (initially detected by Gu04) 
found just above the poles of BLL These are difficult to ob¬ 
serve in the composite HST and [N ii] images but the long-slit 
spectrum of Gu04 at RA. 8° (where the slit goes through the 
symmetry axis and knots) of Fig. shows these knots to be 
present. In fact, the synthetic and observed spectra could not 
be correctly correlated without the inclusion of these compo¬ 
nents. Their kinematical ages, assuming particular parameters, 
indicates that they could have formed around the same time as 
BL2 and EE. Their origination, however, remains unclear. 

BL3, on the other hand, was a difficult component to model 
because of the multitude, direction, and size of the stringy, 
ray-like features that emanate from the nebula’s central region. 
Modelling these rays would not have been possible without the 
use of the detailed spectra from Gu04. What also makes BL3 
interesting is the velocity by which these features expand. The 


rays were modelled with a Hubble-like expansion, which agrees 
very well with the observed spectra. We assumed a radial motion 
and an expanding velocity that is proportional to the distance 
from the nebular centre. As for the density distribution of each 
ray, we performed this visually using the observed spectra and 
HST image in Eig.|^ Gu04 and SG04 described BL3 as probably 
the oldest part of the nebula because of its size, velocity, and its 
somewhat older kinematical age to those of BLl, BL2 and EE. 
Our analysis of BL3 further confirms their findings. 

We initially thought that the equatorial ellipse (EE), which 
is very faint in the composite HST image, might be the earliest 
ejection event due to its large size and halo-like appearance, but 
we obtained a similar kinematical age to Gu04, which puts this 
component in the same age category as BL2 and the knots. This 
suggests that the knots, BL2, and EE could be related to each 
other due to the similarity of their ages. 


4.3. The morpho-kinematics of Hen 2-104. 


In contrast to the results of previous investigations by Corradi 
|et al.| ( |2001[ ) and [Santander-Garcfa et al.| ( |2008| l, we find that 
there is a considerable difference in the kinematical age of the 
two sets of lobes. This is manifest from the two different values 
for the constant of proportionality k that determines the velocity 
field (see Table [^. We find that the ratio between the veloci¬ 
ties for the inner and outer pair of lobes is 6/2.5 = 2.4, and is 
roughly consistent with the ratio of width between them. The ra¬ 
tio between the lengths, as given by the beginning and end of the 
region of knots and rings for the large and small lobes, respec¬ 
tively, is larger. We therefore conclude that the two sets of lobes 
are not coeval as proposed by the earlier studies, and suggest that 
the smaller set of lobes are roughly 2.4 times younger than the 
larger set. 

The value of the constant k - 10 km s * arcsec ' for the red 
knot that we identify with the jet corresponds to a kinematical 
age that is even younger than that of the smaller (inner) lobes. 
However, further study of the outer jets should be performed to 
test for their kinematical age. 

The new detailed reconstruction of the distribution of knots, 
including the regions that are tangential to the line of sight, based 
on the long-slit spectra, suggests that the overall distribution of 
the knots of Hen 2-104 do not appear to be uniform over the 
azimuthal angle of the lobes. There are large regions with con¬ 
siderably fewer knots than in others. There is, however, no clear 
symmetry pattern between the two lobes regarding these varia¬ 
tions in number density of knots and filaments. This result re¬ 
mains tentative, since the spectroscopic coverage in the tangen¬ 
tial regions needs to be complete to make sure that no knots and 
filaments have been missed or misplaced. A better spectroscopic 
coverage with future observations in the tangential regions might 
help to strengthen or weaken this result. 

Erom the modelling of the overall structure of the bipolar 
lobes, in the transition region from the smooth to filamentary 
emission distribution, we found a discontinuity in the opening- 
angle of the lobe. At the position where the knots and filaments 
start, the opening-angle increases suddenly by roughly 10° to 
20° (see the mesh images for an example of this). This change 
in structure hints towards a change in thermal and velocity struc¬ 
ture in the lobes, where a transition from a stable to an unstable 
shock structure happened. We speculate that the low-brightness 
region and brighter rings in the inner bipolar lobes might corre¬ 
spond to similar separate regions. In the future, the brighter rings 
might also break up into smaller scale knots and filaments due 
to thermal and/or dynamic instability. 
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Both on the red and blue side of the lobes we found a conden¬ 
sation that did not fit in the position-velocity profile established 
by the hourglass pattern. Each condensation is close to the sym¬ 
metry axis of the nebula, and are about half-way from the centre 
to the edges of the inner hourglass lobes, approximately along 
the axis. It turns out that assuming they are located along the 
axis with a velocity of 105 km s ' produces the correct place¬ 
ment in the synthetic spectra. Therefore, they might be residuals 
of the collimated jet ejecta. 

Finding the velocity of the outer jet features might also help 
to clarify whether they were ejected in a short burst or in a long 
duration process. In the case of a short burst, the space veloc¬ 
ity of the outer jets should be proportional to their distance and 
about 3-4 times higher than the inner jet condensation. However, 
if a more continuous jet process produces these features, the ve¬ 
locities of the inner knot and outer jet features should be compa¬ 
rable. 


4.4. Separating symbiotic from pianetary nebuiae. 


The main reason for SySts and PNe being distinguished from 
each other in the optical diagram, as proposed by [Gutierrez^ 
Moreno et al. (|1995|l, is their difference in physical conditions 


especially the significantly lower core density in PNe. The same 
is true for all 35 recently discovered SySts in M 31 (see Fig. 5 
from Mikolajewska et al.|2014| l. 

Due to their high densities, young PNe and SySts can oc¬ 
cupy the same region in this diagram, as already pointed out by 
IGutierrez-Moreno et al.| ( [T995| l and [Pereira & Miranda] ( |2005| l. 
Also, the fact that at least two known SySts (Galactic SySt 
Hen 2-147 and extragalactic ICIO SySt-1; Mikolajewska et al. 
1997] Gonsalves et al. 2008 |2015| l are placed in the locus of 


the evolved PNe in the diagram, proves that there are pros and 
cons in the use of such diagnostic diagrams. So we should take 
the results indicated by them with caution. It is however worth 
noting that in the optical spectra of Hen 2-147 and ICIO SySt-1, 
as well as those of M 2-9 and Mz 3, the O vi Raman scattered 
emission line at AA 6830, 7088 were never detected. If so, no 
matter what the other properties of these nebulae are, they would 
be directly classified as SySts; since the presence of this emis¬ 
sion line is a phenomenon uniquely found in SySts, not in all, 
but in at least 50% of them ( Schmid & Schild|1994| ). Altogether 
these arguments seem to indicate that not only two (possibly and 
likely) SySts in NGC 205, but also other known SySts (Galactic 
and extragalactic) are misplaced in the optical diagram of Fig. 
probably due to other (unknown) effects than the density. 

Because of the obvious relevance of the NIR spectrum range 
to study the properties of the red-giant companion in a SySt, 
the way out of the above limitation on the use of purely opti¬ 
cal spectrum is to use it simultaneously with the 2MASS dia¬ 
gram. Though a few PNe appear inside the locus of SySts in 
the 2MASS diagram of Fig. only one SySt is outside the ap¬ 
propriate locus. Also, M 2-9 and Mz 3 are, once again, located 
well inside the locus of SySts, further suggesting their symbiotic 
nature. 

As mentioned above, the presence of the O vi Raman scat¬ 
tered emission line is a powerful criterion to classify a source 
as a SySt (see, for instance, Belczyhski et al.|2000| l. The broad 
emission features at around 6830 and 7088 A are actually O 
VI AA 1032,1038 lines that are Raman-scattered by atomic hy¬ 
drogen, where the scattering H atom initially in the Is state de- 
excites to the 2s state after scattering. A He ii Raman-scattered 
feature, around 4851 A has been reported in the spectrum of 


the symbiotic nova RR Telescopii by jvan Groningen| ( |1993| ), and 
atomic physics tells us that it is also natural to expect that a 
stronger Raman-scattered He ii line will be found around 6545 
A, blueward of [N ii] 6548 A, in the same object, and probably in 
other symbiotic stars ( jNussbaumer et al.|1989 1 . Fee et al. ( [2001 | l 
pointed out that the He ii emission lines require much less H i 
column densities for Raman scattering than the O vi, and that we 
may expect that the former Raman scattering may be observed 
in broader classes of objects than symbiotic stars. 

In the particular case of M 2-9, the He n Raman scattering 
around 6545 A was detected, and the broad Ha wings of this 
nebula was studied and attributed, as much as in known SySts 
(RR Tel and He 2-106) to the hypothesis of having a Raman- 
scattering origin ( Fee et al.|2001| l. As the latter authors discuss, 
at least one young PNe also presents the He n Ram an scattering 


around 4851 A, NGC 7027 (Pequignot et al. 


19971. 


Furthermore, the broad Ha lines with a double or even single 
peak profile ( Van Winckel et ar]|1993| l may be useful in under¬ 
standing the physical conditions at the core of M 2-9, Mz 3, 
Hen 2-104, and even M 1-91. Hen 2-104 shows a broad single 
profile similar to those of known D-type SySts such as RR Tel, 
Rx Pup, and Hen 2-171 (Van Winckel et al. 1993 1 , whereas M 2- 
9 and M 1-91 show double-peaked Ha profiles similar to S-type 
SySts. 

Due to the difference in evolutionary stages between PNe 
and SySts, different morphological characteristics may be ex¬ 
pected, such as in the jets that are observed in both type of sys¬ 
tem. A study of their kinematical ages show that they, more often 
than not, predate the kinematical ages of bipolar lobes associ¬ 
ated with PNe ( Jones et al.|2014| l, such that we would expect 
younger jets to be observable in SySts. This hypothesis is sup¬ 
ported in this work by the detection of a jet feature close to the 
binary centre of Hen 2-104. For these binary cores in SySts it is 
not unusual for long periods > 10 years. For SySts and PNe with 
binary cores, depending on mass ratios and rates of stellar evo¬ 
lution of individual systems, will evolve into nova systems or 
the appropriate analogous sub-type of cataclysmic variable for 
their characteristics. However, further investigation is still nec¬ 
essary to prove these hypotheses and provide a robust technique 
for separating SySts from young PNe. 


5. Conclusions 

We have carried out a morpho-kinematic and spectroscopic 
study of M 2-9, Mz 3, and Hen 2-104 in order to further un¬ 
derstand their structures, kinematics, and nature. The morphol¬ 
ogy and kinematics of these objects have been accurately con¬ 
strained by the use of detailed long-slit spectra. The velocities of 
the individual structures for each object have been determined, 
and a velocity field was fitted to these to reveal the red- and blue- 
shifted motion of their material. Kinematical ages have also been 
obtained for the separate components of M 2-9 and Mz 3. 

The laterally rotating features of M 2-9 were modelled as 
surface brightness enhancements, which are most likely the re¬ 
sult of excitation on the walls of the bulb component by a col¬ 
limated spray of high velocity particles. On the other hand, the 
morpho-kinematic modelling of Mz 3 has revealed a second in¬ 
ner bipolar structure within BFl and a set of knots found just 
above the apexes of BFl. For Hen 2-104, we found what ap¬ 
pears to be two knotty jet structures close to the symmetry axis 
and found at ~11" and 13" on either side of the outer bipolar 
lobes. Also, according to our modelling of Hen 2-104, the kine- 
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matical age of the inner and outer bipolar lobes are different, 
with the outer lobes 2.4 times older. 

Apart from just studying the morphological and kinematical 
characteristics, shape also provided us with the density distribu¬ 
tion of each nebula in the form of an ASCII data cube, which 
will allow for the future work of creating 3-D photo-ionisation 
models. These data cubes can be used as an input for other ex¬ 
ternal codes such as moccasin ( |Ercolano et al.||2b03| l to create 
complex 3-D photo-ionisation models and investigate, in detail, 
the chemical structure of the nebulae, while avoiding the time 
consuming task of reproducing the density distribution of these 
models from the beginning. 

The optical emission line ratios and 2MASS colour values 
present further evidence for the symbiotic nature of M 2-9 and 
Mz 3. Both diagnostic diagrams have been shown to success¬ 
fully separate a multitude of known SySts from genuine PNe. 
However, the NIR colour-colour diagram can not be used for the 
classihcation of SySts since several other objects such as Be, 
B[e], WR, young PNe, and H ii regions, share the same range of 
colours. 

The Her prohle of M 2-9 shows a broad line with an absorp¬ 
tion dip and two asymmetric components that are usually found 
in SySts because of the fast stellar wind and the high mass-loss 
rate of the cold companion. This is indicative of a binary symbi¬ 
otic system. 
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Appendix A: Tables and Mesh Images 

The optical emission line intensities for the various SySts and 
PNe are shown in Tables. |^and|^respectively. The 2MASS NIR 
colour magnitude values for the same sample of SySts and PNe 
are presented in Tables. and |7] Shown in Fig. [^are the 3-D 
SHAPE images of M 2-9, Mz 3, and Hen 2-104 respectively. In all 
of these images, the red line represents the top view, the blue line 
represents the right (or side) view, and the green line represents 
the symmetry (or major) axis. 
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Table 4. Line intensities for the nucleus of specific symbiotic nebulae from various literature 


Object Type [Hr]4340 A [O iii]4363 A [H/3]4861 A [O iii]5007 A 


AS210 

SySt 

64.3 

Hen 2-38 

SySt 

53.9 

RRTel 

SySt 

50.6 

RXPup 

SySt 

54.5 

HMSge 

SySt 

45.2 

He 2-106 

SySt 

51.2 

He 2-127 

SySt 

45.0 

He 2-171* 

SySt 

34.3 

H 1-36* 

SySt 

46.1 

SS 122 

SySt 

56.9 

He 2-25 

SySt 

44.4 

Th2-B 

SySt 

36.2 

He 2-171 

SySt 

45.1 

He 2-127 

SySt 

46.3 

SS 122 

SySt 

35.0 

He 2-106 

SySt 

44.8 

Hen 2-104 ^ 

SySt 

39.7 

Hen 2-104 ^ 

SySt 

53.9 

He 2-171 

SySt 

45.0 

SS 122 

SySt 

38.6 

H2-38 

SySt 

47.3 

H2-5 

SySt 

41.2 

H 1-36** 

SySt 

45.1 

H 1-25 

SySt 

40.2 

Hen 3-1342* 

SySt 

44.6 

H2-34 

SySt 

36.8 

Hen 3-1591 

SySt 

49.9 

Hen 3-1761 

SySt 

54.2 

Hen 3-1410 

SySt 

43.1 

Th 3-29 

SySt 

45.0 

Wray 16-377 

SySt 

44.1 

SS 96 

SySt 

38.9 

SS 71 

SySt 

41.3 

AS210 

SySt 

44.5 

AS269 

SySt 

37.2 

RRTel 

SySt 

48.3 

K66 

SySt 

30.8 


77.0 

100 

171.1 

220.1 

100 

707.7 

42.8 

100 

193.8 

82.2 

100 

234.4 

71.4 

100 

430.9 

43.1 

100 

94.3 

32.8 

100 

115.4 

25.8 

100 

384.8 

55.4 

100 

1453.0 

47.4 

100 

154.3 

161.6 

100 

475.9 

28.0 

100 

104.4 

95.0 

100 

446.6 

50.8 

100 

214.7 

30.8 

100 

150.9 

53.9 

100 

308.5 

51.0 

100 

187.5 

71.8 

100 

301.9 

44.6 

100 

336.4 

23.0 

100 

117.2 

113.2 

100 

417.5 

24.1 

100 

34.9 

52.0 

100 

1492.7 

21.1 

100 

43.9 

1.4 

100 

11.0 

20.0 

100 

52.8 

118.7 

100 

281.2 

35.7 

100 

57.1 

12.9 

100 

40.5 

28.8 

100 

39.6 

14.7 

100 

20.6 

17.3 

100 

43.9 

11.9 

100 

36.6 

16.8 

100 

93.6 

56.5 

100 

7.7 

33.9 

100 

155.2 

14.3 

100 

61.2 


Reference 


de Freitas Pacheco & Costa (1992 1 


de Lreitas Pacheco & (Josta 


de Lrettas Pacheco & (Josta 


de Freitas Pacheco & (Josta 


de Freitas Pacheco & (Josta 


(Josta & de Lreitas Pacheco 


(Josta & de Lreitas Pacheco 


(Josta & de Lrettas Facheco 


(Josta & de Lreitas Facheco 


(Josta & de Lreitas Pacheco' (19941 


(T952i 


iTWZi 


(19921 


(T992i 


(T995i 


(T995i 


(199211 


(199211 



(Jorradi 

(1995 

1 


(Jorradi 

(1995 

1 


Pereira 

1995 


Pereira 

1995 

Pereira 

1995 

Pereira 

1995 

Pereira 

T995 



Luna & (Josta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 

Luna & Losta 

(2005 


Notes. A list of symbiotics chosen from literature with extinction correction intensity, ^objects of interest within this paper, *objects bordering 
between regions in the optical diagnostic diagram, **objects that fall outside their associated (SySt) region. 
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Table 5. Line intensities for the nucleus of specific PNe from various literature 


Object 

Type 

[H7]4340 a 

[O iii]4363 A 

[1^14861 A 

[O Iii]5007 A 





Reference 





MyCn 18 

PN 

46.8 

0.45 

100 

32T 


King 

sburgh & Barlow 

(1994i 

Mz 3^* 

PN 

52.5 

3.984 

100 

32.58 




Zhang & Liu 

(300 

3i 



Mz3l' 

PN 

46.2 

4.963 

100 

44.96 





Smith 

30 





NGC 1535 

PN 

46.8 

12.65 

100 

1189.38 



Krabbe & Copetti 

2006 

1 


NGC 2438 

PN 

46.8 

9.56 

100 

1045.42 



Krabbe & Copetti 

2006 

1 


NGC 2440 

PN 

46.8 

23.98 

100 

1312.31 



Krabbe & Copetti 

2006 

1 


NGC 3132 

PN 

46.8 

4.56 

100 

885.72 



Krabbe & Copetti 

2006 

1 


NGC 3242 

PN 

46.8 

14.05 

100 

1204.84 



Krabbe & Copetti 

2006 

1 


NGC 6302 

PN 

46.8 

36.46 

100 

1271.16 



Krabbe & Copetti 

2006 

1 


NGC 7009 

PN 

46.8 

7.57 

100 

1170.75 



Krabbe & Copetti 

2006 

1 


NGC 7027 

PN 

47.6 

25.4 

100 

1657.9 




Zhang et al. 

(pD 

07 1 



M 2-9 '' 

PN 

53.1 

46.8 

100 

268.4 


Torres-Peimbert et al. 

(2010 

1 

M 1-91 

PN 

45.9 

36.1 

100 

232.5 


Ibrres-Peimbert et al. 

(2010 

1 

NGC 2392 

PN 

39.35 

14.9 

100 

903.5 




Zhang et al. 

pUT 

L \ 



PB5 

PN 

47.17 

22.39 

100 

1422.24 





Gorny 


)T^ 

1 




PB 3 

PN 

48.24 

14.26 

100 

1306.84 





Gorny 

(THET 

1 




He 2-36 

PN 

47.26 

19.46 

100 

1345.42 





Gorny 

(2U14 

1 




He 2-39 

PN 

47.24 

22.88 

100 

1314.91 





Gorny 

(2U14 

1 




He 2-63 

PN 

47.09 

12.67 

100 

1020.17 





Gorny 

(2U14 

1 




He 2-67 

PN 

47.28 

6.55 

100 

1165.17 





Gorny 

(2U14 

1 




He 2-70 

PN 

47.15 

34.84 

100 

1022.69 





Gorny 

(2U14 

1 




He 2-84 

PN 

47.95 

11.19 

100 

1121.92 





Gorny 

(2U14 

1 




He 2-85 

PN 

47.03 

14.81 

100 

1365.70 





Gorny 

(2U14 

1 




He 2-86 

PN 

46.64 

3.42 

100 

843.42 





Gorny 

(2U14 

1 




He 2-96 

PN 

46.60 

3.08 

100 

843.37 





Gorny 

(2U14 

1 




My 60 

PN 

47.90 

19.46 

100 

1196.34 





Gorny 

(2U14 

1 




Pe 1-1 

PN 

48.26 

7.67 

100 

944.27 





Gorny 

(2U14 

1 




ESO 095-1 

PN 

47.13 

20.88 

100 

1538.56 





Gorny 

(2U14 

1 




ESO 320-2 

PN 

47.46 

9.71 

100 

349.44 





Gorny 

(2U14 

1 




BMPJl 128-6 

PN 

47.32 

37.08 

100 

1352.23 





Gorny 

(Sura 





NGC 3918 

PN 

46.65 

20.98 

100 

1578.26 





Gorny 

(2U14 

1 




K 1-23 

PN 

46.96 

8.48 

100 

934.96 





Gorny 

(2U14 

1 




Wray 16-128 

PN 

47.63 

11.32 

100 

440.80 





Gorny 

(2BT7 





Me We 1-3 

PN 

47.46 

18.61 

100 

676.26 





Gorny 

(2U14 

1 




MaC 1-2 

PN 

47.03 

15.85 

100 

1490.70 





Gorny 

(2U14 

1 




He 2-102 

PN 

46.74 

8.33 

100 

1000.89 





Gorny 

(Sura 





He 2-103 

PN 

52.47 

9.72 

100 

802.13 





Gorny 

(2U14 

1 




He 2-109 

PN 

46.93 

9.04 

100 

1087.32 





Gorny 

(2U14 

1 




He 2-115 

PN 

46.68 

2.58 

100 

558.55 





Gorny 

(2BK 





He 2-117 

PN 

46.54 

2.22 

100 

721.78 





Gorny 

(2U14 

1 




He 2-133 

PN 

52.81 

8.08 

100 

1183.33 





Gorny 

(2U14 

1 




NGC 5939 

PN 

47.27 

17.28 

100 

946.07 





Gorny 

(217T4 





NGC 6326 

PN 

46.52 

14.14 

100 

1330.61 





Gorny 

(2U14 

1 




StWr4-10 

PN 

47.13 

10.95 

100 

814.49 





Gorny 

(2U14 

1 




K3-31 

PN 

50.25 

13.13 

100 

1466.56 

Garcia-He 

rnande 2 

t&G( 

3rny 

( 

2014 

K3-56 

PN 

47.38 

11.71 

100 

506.08 

Garcia-Hernandez & Gorny 

( 

2014 

K3-87 

PN 

47.30 

11.11 

100 

558.00 

Garcia-Hernandez & Gorny 

( 

2014 

K4-41 

PN 

46.65 

5.63 

100 

994.52 

Garcia-Hernandez & Gorny 

( 

2014 

M2-50 

PN 

47.09 

14.67 

100 

1190.65 

Garcia-Hernandez & Gorny 

( 

2014 

B1 2-1 

PN 

52.03 

6.44 

100 

981.25 

Garcia-Hernandez & Gorny 

( 

2014 

Hen 2-440 

PN 

46.08 

1.96 

100 

437.94 

Garcia-Hernandez & Gorny 

( 

2014 

NGC 6807 

PN 

46.86 

13.19 

100 

1353.14 

Garcia-Hernandez & Gorny 

( 

2014 

M 1-69 

PN 

46.78 

7.12 

100 

1216.84 

Garcia-Hernandez & Gorny 

( 

2014 

IC 4846 

PN 

46.84 

7.31 

100 

1046.40 

Garcia-Hernandez & Gorny 

( 

2014 

K4-8 

PN 

46.90 

10.15 

100 

1235.16 

Garcia-Hernandez & Gorny 

( 

2014 

Sa 2-237 

PN 

46.99 

10.47 

100 

1070.36 

Garcia-Hernandez & Gorny 

( 

2014 

MaC 1-11 

PN 

49.07 

7.64 

100 

869.54 

Garcia-Hernandez & Gorny 

( 

2014 


Notes. A list of planetary nebulae (PN) chosen from literature with extinction correction intensity, ^objects of interest within this paper, *objects 
bordering between regions in the optical diagnostic diagram. “ the intensity of the [O iii]5007 A line is obtained from the [O iii]4959 A line 
multiplied by 3. 
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Table 6. 2MASS NIR (HJKs) colour magnitude values of various SySts. 


Object 

Type 

H 

J 

Ks 

He 2-171 

SySt 

9.679 

8.028 

6.640 

He 2-127 

SySt 

10.585 

9.433 

8.530 

SS 122** 

SySt 

9.958 

9.796 

9.567 

Hen 2-104 '' 

SySt 

11.005 

8.957 

7.145 

He 2-106 

SySt 

9.584 

7.756 

6.135 

He 2-171 

SySt 

9.679 

8.028 

6.640 

H2-38 

SySt 

8.266 

7.240 

6.478 

H2-5 

SySt 

7.161 

6.030 

5.556 

H 1-36* 

SySt 

11.771“ 

9.873 

7.625 

H 1-25 

SySt 

10.415 

8.633 

6.949 

Hen 3-1342 

SySt 

9.605 

8.693 

8.348 

H2-34 

SySt 

9.423 

8.125 

7.607 

Hen 3-1591 

SySt 

10.504 

9.680 

9.073 

Hen 3-1761 

SySt 

6.681 

5.797 

5.446 

Hen 3-1410 

SySt 

10.315 

9.158 

8.572 

Th 3-29 

SySt 

8.924 

7.572 

6.921 

Wray 16-377 

SySt 

9.802 

8.914 

8.605 

SS73 96 

SySt 

8.187 

6.905 

6.421 

SS73 71 

SySt 

11.374 

10.569 

10.162 

AS210 

SySt 

9.087 

7.555 

6.256 

AS269 

SySt 

12.148 

11.308 

9.936 

RRTel 

SySt 

7.302 

6.079 

4.902 

K6-6 

SySt 

9.423 

8.199 

7.523 

He 2-25 

SySt 

13.412 

12.115 

10.406 

Th 2-B 

SySt 

12.951 

11.405 

9.891 

RXPup 

SySt 

6.397 

4.981 

3.782 

HM Sge 

SySt 

7.620 

6.020 

4.654 


Notes. The data were obtained from the 2MASS All-Sky Catalog of Point Sources i Cutri et al. 2003b i. " uncertain value, objects of interest within 
this paper, *objects bordering between regions in the 2MASS colour-colour diagram, **objects that fall outside their associated (SySt) region. 
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Table 7. 2MASS NIR (HJKs) colors of various PNe. 


Object 

Type 

H 

J 

Ks 

NGC 1535 

PN 

12.539 

12.647 

12.577 

NGC 2438 

PN 

12.672 

12.526 

12.509 

NGC 2440 

PN 

15.671 

15.037 

15.203 

NGC 3132 

PN 

9.754 

9.731 

9.723 

NGC 3242 

PN 

12.102 

12.461 

11.748 

NGC 6302 

PN 

11.255 

11.706 

9.442 

NGC 7009 ' 

PN 

9.325 

9.456 

8.748 

NGC 7027 

PN 

9.613 

8.37 

7.443 

NGC 2392 

PN 

10.872 

10.919 

10.936 

MyCn 18 

PN 

11.554 

11.603 

10.808 

Mz3'' 

PN 

9.352 

7.355 

5.608 

PB5 

PN 

12.488 

11.897 

10.625 

PB 3 

PN 

13.899 

13.703 

12.651 

He 2-36 

PN 

9.867 

9.623 

9.463 

He 2-39 

PN 

13.474 

12.614 

12.338 

He 2-63 

PN 

15.55 

15.427 

14.613 

He 2-67 

PN 

13.43 

13.431 

12.454 

He 2-70 

PN 

12.544 

11.497 

11.180 

He 2-84 

PN 

14.527 

13.422 

13.077 

He 2-85 

PN 

15.198 

14.305 

13.681 

He 2-86 

PN 

12.155 

11.853 

10.789 

He 2-96 

PN 

12.577 

12.228 

11.178 

My 60" 

PN 

12.604 

12.379 

11.635 

Pe 1-1 

PN 

12.457 

12.065 

10.931 

ESO 095-1 

PN 

15.201 

13.798 

11.867 

ESO 320-2 

PN 

15.166 

14.157 

13.606 

BMP J1128-6121 

PN 

14.155 

13.478 

12.340 

NGC 3918 

PN 

11.691 

11.568 

10.422 

K 1-23 

PN 

16.701 

15.941 

16.584 

Wray 16-128 

PN 

15.761 

15.078 

14.872 

MaC 1-2 

PN 

13.713 

13.446 

13.315 

He 2-102 

PN 

14.636 

14.31 

13.661 

He 2-103 

PN 

14.971 

14.34 

13.612 

He 2-108 

PN 

12.173 

12.103 

11.877 

He 2-109 

PN 

16.338 

15.832 

15.061“ 

He 2-115 

PN 

12.033 

11.724 

10.598 

He 2-117 

PN 

11.811'* 

11.424“ 

10.195“ 

He 2-133 

PN 

12.132 

11.72 

10.579 

NGC 5939 

PN 

14.777 

13.782 

12.700“ 

NGC 6326 

PN 

14.747 

14.939 

14.104 

StWr4-10 

PN 

14.335 

14.209 

13.824 

M 2-9'' 

PN 

11.198 

9.177 

6.996 

M 1-91 

PN 

14.099 

11.841 

9.710 

K3-31 

PN 

13.806 

13.339 

12.300 

K3-56 

PN 

15.973 

14.895“ 

14.689 

K3-87 

PN 

16.214 

15.684 

14.919 

K4-41 

PN 

15.602** 

13.796“ 

14.473 

M2-50 

PN 

15.333 

15.242 

14.456 

B1 2-1 

PN 

13.473 

12.607 

11.16 

Hen 2-440 

PN 

13.603 

12.686 

12.414 

NGC 6807 

PN 

12.753 

12.826 

12.069 

M 1-69 

PN 

13.261 

13.207 

12.286 

IC 4846 

PN 

12.681 

12.835 

11.897 

K4-8 

PN 

14.792 

14.802 

13.946 

Sa 2-237 

PN 

13.065 

12.566 

11.817 

MaC 1-11 

PN 

14.179 

13.779 

13.107 


Notes. The data have been obtained from the 2MASS All-Sky Catalog of Point Sources ^|Cutri et ^ 2003b| . ^ objects of interest within this paper. 
^"from[Skrutskie et al.|(|200^. “ uncertain value. 
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Hen 2-104 
Front view 


Fig. 9. Mesh images of M 2-9, Mz 3, and Hen 2-104 showing their components and sub-components. Shown on the left column 


are the render views for the images seen in Sects. [3T||3.2| and 3.3 
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